One-dimensional superfiuid Bose-Fermi mixture: mixing, demixing and bright solitons 
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We study a ultra-cold and dilute superfluid Bose-Fermi mixture confined in a strictly one- 
dimensional atomic waveguide by using a set of coupled nonlinear mean-field equations obtained 
from the Lieb-Liniger energy density for bosons and the Gaudin-Yang energy density for fermions. 
We consider a finite Bose-Fermi inter-atomic strength g^f and both periodic and open boundary 
conditions. We find that with periodic boundary conditions, i.e. in a quasi-lD ring, a uniform 
Bose-Fermi mixture is stable only with a large fermionic density. We predict that at small fermionic 
densities the ground state of the system displays demixing if gbj > and may become a localized 
Bose-Fermi bright soliton for gi,f < 0. Finally, we show, using variational and numerical solution 
of the mean-field equations, that with open boundary conditions, i.e. in a quasi-lD cylinder, the 
Bose-Fermi bright soliton is the unique ground state of the system with a finite number of parti- 
cles, which could exhibit a partial mixing-demixing transition. In this case the bright solitons are 
demonstrated to be dynamically stable. The experimental realization of these Bose-Fermi bright 
solitons seems possible with present setups. 
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I. INTRODUCTION 



Effects of quantum statistics are strongly enhanced 
in strictly one-dimensional (ID) systems These ef- 
fects can be investigated with ultra-cold and dilute gases 
of alkali-metal atoms, which are now actively studied 
in the regime of deep Bose and Fermi degeneracy 0. 
Recently two experimental groups 0, Ql have reported 
the observation of the crossover between a ID quasi 
Bose-Einstein condensate (BEC) in the weak-coupling 
mean-field Gross-Pitaevskii (GP) domain and a Tonks- 
Girardeau (TG) gas [1, 0] with ultra-cold ®^Rb atoms 
in highly elongated traps. A rigorous theoretical analy- 
sis of the ground-state properties of a uniform ID Bose 
gas was performed by Lieb and Liniger (LL) forty four 
years ago An extension of the LL theory for finite 
and inhomogeneous ID Bose gases has been proposed 
on the basis the local density approximation (LDA) 
The LDA is improved by including a gradient term that 
represents additional kinetic e nerg y associated with the 
inhomogencity of the gas [1, i, [la lill, E, [ll- 

In the last few years several experimental groups 
have observed the crossover from the Bardeen-Cooper- 
Schricffer (BCS) state of Cooper Fermi pairs to the 
BEC of molecular dimers with ultra-cold twq-hypcrfinc- 
component Fermi vapors of atoms [l^jSjIill and ^Li 
atoms [13, [111 ■ It is well known that purely attractive po- 
tentials have bound states in ID and 2D for any strength, 
contrary to the 3D case [l^ . A rigorous theoretical anal- 
ysis of the ground-state properties of a uniform ID Fermi 
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gas was performed by Gaudin and Yang (GY) forty years 
ago [lOl ■ For repulsive interaction the GY model gives a 
Tomonaga-Luttinger liquid [2]| . while for attractive in- 
teraction it describes a Luther-Emery (superfluid) liquid 
[2^ . The ground state of a weakly attractive ID Fermi 
gas is a BCS-like state with Cooper pairs, whose size 
is much larger than the average inter-particle spacing. 
The strong coupling regime with tightly bound dimers 
is reached by increasing the magnitude of the attractive 
inter-atomic strength. In this regime the fermion pairs 
behave like a hard core Bose gas (TG gas), or equiva- 
lently like ID noninteracting fermions [23l. |24|. 

Degenerate Bose-Fermi mixtures of alkali-metal atoms 
have been experimentally observed in ^^^Li 0, ®Li- 
23Na [p, and ^OR-S^Rb [2|. In these mixtureg^ the 
theoretical investigation of phase separation [2^, [30, [Ml, 
m, [H, [13, [11] and soliton-like structures has drawn sig- 
nificant attention. Bright solitons have been observed in 
BEC's of Li [111 and Rb [13] atoms and studied subse- 
quently nil . It has been demonstrated using microscopic 
[3^ and mean-fleld hydrodynamic models that the 
formation of stable fermionic bright and dark solitons is 
possible in a degenerate Bose-Fermi mixture as well as 
in a Fermi-Fermi mixture [JTI with the fermions in the 
normal state in the presence of a sufficiently attractive 
interspecies interaction which can overcome the Pauli re- 
pulsion among fermions. The formation of a soliton in 
these cases is related to the fact that the system can 
lower its energy by forming high density regions (soli- 
tons) when the interspecies attraction is large enough to 
overcome the Pauli repulsion in the degenerate Fermi gas 
(and any possible repulsion in the BEC) [i^l- There have 
also been studies of mixing-demixing transition in degen- 
erate Bose-Fermi [43| and Fermi- Fermi [4^ mixtures with 
fermions in the normal state. 

After observing the degenerate Fermi gas and the 
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realization of BCS condensed superfluid phase of the 
fermionic system [l^, the BCS-Bose crossover [4^ in it 
seems to be under control [13, E, [H, [13, [H] by ma- 
nipulating the Fermi-Fermi interaction near a Feshbach 
resonance [47t by varying a uniform external magnetic 
field. Naturally, the question of the possibility of bright 
solitons in a superfluid fermionic condensate assisted by 
a EEC (with attractive interspecies Bosc-Fcrmi interac- 
tion) needs to be revisited. Although the fermionic BCS 
or molecular dimer phase is possible with an attractive 
Fermi-Fermi attraction, once formed such a phase is in- 
herently repulsive and self-defocussing and will not lead 
to a bright soliton spontaneously. 

In this paper we investigate a ID superfluid Bosc-Fcrmi 
mixture composed of bosonic atoms, well described by 
the LL theory, and superfluid fermionic atoms, well de- 
scribed by the GY theory with the special intention of 
studying the localized structures or bright solitons in this 
superfluid mixture. We derive a set of coupled nonlinear 
mean-field equations for the Bose-Fermi mixture which 
we use in the present study. The solution of this mean- 
field equation is considered two types of boundary con- 
ditions: periodic and open. It has been recently shown 
jisl . [49I . [soj that an attractive BEC with periodic bound- 
ary conditions, which can be experimentally produced 
with a quasi-lD ring [5l|, displays a quantum phase 
transition from a uniform state to a symmetry-breaking 
state characterized by a localized bright-soliton conden- 
sate (soj . Here we show that a similar phenomenon 
appears in the ID Bose-Fermi mixture with an attrac- 
tive Bose-Fermi scattering length. Instead, with a re- 
pulsive Bose-Fermi scattering length or with the increase 
of the number of Fermi atoms leading to large repulsion 
we find a phase separation between bosons and fcrmions 
(a mixing-demixing transition), in analogy with previous 
theoretical [1^, HO, HU, i32„ j53„ ^ and experimental 
studies [H with 3D trapped ^Li-'^Li and ^^K-'^°K mix- 
tures. Finally, we predict that with open boundary condi- 
tions, i.e. in a inflnite quasi-lD cylinder, the ground-state 
of the mixture with sufficiently attractive Bose-Fermi in- 
teraction is always a localized bright soliton. 

The paper is organized as follows. In Sec. II we de- 
scribe the ID model we use in our investigation of the su- 
perfluid Bose-Fermi mixture. The Lagrangian density for 
bosons is appropriate for a TG to BEC crossover through 
the use of a quasi-analytic LL function. The fermions are 
treated using the GY model through a quasi-analytic GY 
function. The Bose-Fermi interaction is taken to be a 
standard contact interaction. The Eulcr-Lagrange equa- 
tions for this system arc a set of coupled nonlinear mean- 
field equations which wc use in the present study. In Sec. 
Ill we consider the system with periodic boundary con- 
ditions and obtain the thresholds for demixing, for the 
formation of localized bright solitons (for a sufficiently 
strong Bose-Fermi attraction), and for the existence of 
mixing, i.e. states with constant density in space (for 
weak Bose-Fermi attraction and for Bosc-Fcrmi repul- 
sion). We also present a modulational instability analy- 



sis of a uniform solution of the mixture and obtain the 
condition for the appearance of bright solitons by mod- 
ulational instability. In Sec. IV we consider the mixture 
with open boundary conditions and derive a variational 
approximation for the solution of the mean-field equa- 
tions of the model using a Gaussian- type ansatz. We 
investigate numerically the bright Bosc-Fcrmi solitons, 
demonstrate their stability under perturbation, and com- 
pare the numerical results with those of the variational 
approach. Finally, in Sec. V we present a summary and 
discussion of our study. 



II. THE MODEL 

We consider a mixture of iV;, atomic bosons of mass ni}, 
and Nf superfluid atomic fcrmions of mass m/ at zero 
temperature trapped by a tight cylindrically symmetric 
harmonic potential of frequency uj± in the transverse (ra- 
dial cylindric) direction. We assume factorization of the 
transverse degrees of freedom. This is justified in ID 
confinement where, regardless of the longitudinal behav- 
ior or statistics, the transverse spatial profile is that of 
the single-particle ground-state [13, H^l . The transverse 
width of the atom distribution is given by the character- 
istic harmonic length of the single-particle ground-state: 
a±j = ^yh/{^mjU!±), with j = b, f. The atoms have an ef- 
fective ID behavior at zero temperature if their chemical 
potentials are much smaller than the transverse energy 
huj± [33 . [5^ . The boson-boson interaction is described 
by a contact pseudo potential with repulsive (positive) 
scattering length at- The fermion-fermion interaction is 
modelled by a contact pseudo potential with attractive 
(negative) scattering length af. The boson- fcrmion in- 
teraction is instead characterized by a contact pseudo 
potential with scattering length abf, which can be re- 
pulsive or attractive. To avoid the confinement-induced 
53[ we take at, |a/|, \abf \ <C a±b, a±f. 
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We use a mean-field effective Lagrangian to study the 
static and collective properties of the ID Bose-Fermi mix- 
ture. The Lagrangian density C of the mixture reads 



C — Cb + Cf + Cbf ■ 
The term Cb is the bosonic Lagrangian, defined as 
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Cb = ihipl^ijjb 



d_ 



2^b'^''d^'^'' 



2mb 



rribgb 



(2) 



where tpb{z,t) is the hydrodynamic field of the Bose 
gas along the longitudinal axis, such that nb{z,t) = 
|?/'h(z,t)p is the ID local probability density of bosonic 
atoms and Qb = 2huj± Ub is the ID boson-boson interac- 
tion strength (g^ > 0). G{x) is the Licb-Liniger function 
[5^ . defined as the solution of a Frcdholm equation for 



3 



X > 0, and such that 
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for < a; < 1, 
for a; ^ 1. 



(3) 



In the extreme weak-couphng hmit x — > +0, and G{x) — > 
X and consequently, the bosonic Lagrangian above re- 
duces to the standard mean-field GP Lagrangian. In the 
static case the Lagrangian density Cb reduces exactly 
to the energy functional recently introduced by Licb, 
Seiringer and Yngvason [§]. In addition, Cb has been 
successfully used to determine the collective oscillation 
of the ID Bose gas with longitudinal harmonic confine- 
ment In the strong coupling limit x +oo, and 
G{x) ^ 7r^/3 while the Lagrangian above reduces to the 
strongly repulsive bosonic Lagrangian in the TG limit 
@. As' X changes from -1-0 to 4-00 the above Lagrangian 
shows a continuous transition from the GP BEG to TG 
phase. 

The fermionic Lagrangian density Cf is given instead 

by 
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where tpf{z, t) is the field of the ID superfluid Fermi gas, 
such that nf{z,t) = is the ID fermionic local 

density along the longitudinal axis and gf — 2hLj±^ ct/ is 
the ID fcrmion-fcrmion interaction strength (gj < 0). 
F{x) is the Gaudin-Yang function [5^, defined as the 
solution of a Fredholm equation for a; < 0, and such that 

I23,l2i 
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for x <C —1, 
for - 1 < .1 < 0. 



(5) 

In the static and uniform case the Lagrangian density C / 
reduces exactly to Gaudin-Yang energy functional [20| . 
that has been recently used by Fuchs et al. [23| . 

The physical content of the fermionic Lagrangian ^ 
is easy to understand. For x —0 we are in the do- 
main of weak Fermi- Fermi attraction (af,gf cor- 
responding to the BGS hmit), while F{x) 77^/12 
from Eq. Gonsequently, the fermionic interaction 

term in Eq. ([4]) involving the F(x) function reduces to 
(;iV2TO/)(7rVl2)|V'/l'^ (independent of Fermi-Fer mi in- 
teraction strength g/ or Fermi- Fermi scattering length 
a/), which is the Pauli repulsive term of nonintcract- 
ing fermions considered in Rcf. [s^ . [56j arising due to 
the Pauli exclusion principle and not due to the funda- 
mental Fermi-Fcrmi interaction implicit in the scattering 
length ttf via gf. This is quite expected as by taking the 
X —^ —0 limit we switch off the Fermi-Fermi attraction 
and pass from the BGS domain to a degenerate noninter- 
acting Fermi gas (in non-superfluid phase) studied in Ref. 
[5^ corresponding to a; — > -1-0. [Although the probability 



densities in the x — > -1-0 (non-superfluid degenerate Fermi 
gas) and a; ^ — (superfluid BGS condensate) hmits 
are identical, they correspond to different physical states. 
The superfluid phase gas a "gap" associated with it de- 
scribable by the BGS equation.] As x passes from —0 to 
— cx), the Fermi- Fermi attraction increases and we move 
from the superfluid BGS regime of weakly-bound Gooper 
pairs to the unitarity regime of strongly-bound molecular 
fermions. (The unitarity limit corresponds to infinitely 
large Fermi-Fermi attraction: a/ — > —00.) As Fermi- 
Fermi scattering length |a/| increases from the BGS to 
the unitarity limit, the interaction energy in Eq. Q gets 
a term dependent on ay. However, in the unitarity limit 
a/ ^ — 00 and all the fermions will be paired to form 
strongly bound nonintcracting molecules and the Fermi- 
Fermi interaction term in Eq. ^ behaves like that of 
a Tonks- Girardeau gas of bosons — bosonic molecules — 
and again becomes independent of ay [23j . The system 
then becomes a Bose gas of molecules and its interaction 
energy is no longer a function of ay . 

Finally, the Lagrangian density Cbf of the Bose-Fermi 
interaction taken to be of the following standard zero- 
range form [4^, [s^l 



Cbf = -gbf li/'fePlV'/P 



(6) 



where gbf ~ 2TiLo^abf is the ID boson-fermion interaction 
strength. 

The Euler-Lagrangc equations of the Lagrangian C are 
the two following coupled partial differential equations: 
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where rij = ,j = b,f are probability densities for 
bosons and fermions, respectively, with the normalization 
I-^njdz = Nj. 

It is convenient to work in terms of dimensionless 
variables defined in terms of a frequency uj and length 
I = yjhj [rtibiS) by ~ ■'Ai/v^i ^ = 2t/tj, z = zl, 
9j = 5j^^/(2mf,Z), and gbf = gbfh^ /{2mbl). In terms 
of these new variables Eqs. ([7]) and ^ can be written as 
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where we have dropped the hats over the variables, and 
where A = mi,/mf, rij = \ipj\'^,j = b, f are probabihty 
densities for bosons and fcrmions, respectively, with the 
normalization J_^njdz = Nj. However, in the present 



study wc set A = 1 in the following. All results reported 
in this paper are for this case. This should approximate 
well the Bosc-Fcrmi mixtures ''Li-^Li and '^^K-^''K of ex- 
perimental interest. 

The bosonic and fermionic nonlinearities in Eqs. ([9|) 
and (|10|). respectively, have complex structures in gen- 
eral. However, in the weak-coupling GP limit [gt, 0) 
the bosonic nonlinearity in Eq. ([9]) is cubic and turns 
quintic in the strong-coupling TG limit (<?;, — s- +oo). The 
fermionic nonlinearity in Eq. pTH) becomes quintic in 
both weak {gp —0) and strong {gf — > — oo) coupling 
but with coefficients 7r^/12 and 7r^/48, respectively. 

For stationary states the solution of Eqs. ^ and pU]) 
have the form ipj = (pj exp(—ifijt) where fij are the re- 
spective chemical potentials. Consequently, these equa- 
tions reduce to 
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(11) 



(12) 



A repulsive Bosc-Bose interaction is produced by a posi- 
tive gt, while an attractive Fermi- Fermi interaction cor- 
responds to a negative gf. 



III. MIXTURE WITH PERIODIC BOUNDARY 
CONDITIONS 
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FIG. 1: (Color online) Region of stability of the homoge- 
neous mixture denoted by shaded (gray) area in the plane 
(ubjUf). Different panels correspond to different values of 
the Bose- Fermi interaction strength g^f. We set gt — 0.1 and 
(a) gf = —0.1 (corresponding to a weak Fermi- Fermi attrac- 
tion in the BCS phase) and (b) <?/ = —30 (corresponding to a 
strong Fermi-Fermi attraction, with fermions in the unitarity 
regime) . 



A. General Considerations 

Here we consider the ID Bosc-Fermi mixture with pe- 
riodic boundary conditions. These boundary conditions 
can be used to model a quasi-lD ring of radius R. If 
the radius is much larger than the transverse width, i.e. 
R ^ a±i, J = b, /, then effects of curvature can be ne- 
glected [50| and one can safely use the Lagrangian dei isity 
([T]) with z = R9, where 9 the azimuthal angle [sO, l5ll |. 
The energy density of the uniform mixture is immedi- 
ately derived from the Lagrangian density ([T|), with 
(|4]) and ([6]), dropping out space-time derivatives. It is 
given by 

In this case one can have a uniform mixture or the for- 
mation of soliton-like structures and in the following we 



study the condition for these possibilities to happen. Ob- 
viously, for the ID uniform mixture in the ring of radius 

R, we have Ub = Nb/{2T:R) and uf = Nf/{2T:R). 

The uniformly mixed phase is energetically stable if 
its energy is a minimum with respect to small variations 
of the densities n/ and rib, while the total number of 
fermions and bosons are held fixed. To get the equilib- 
rium densities one must then minimize the function 



£{nb,nf) = £{nb,nf) - ij.b{nl,n*f)nb - ^J,f{nl,n*f)nf 

(14) 

where ^b and /i / are Lagrange multipliers (imposing that 
the numbers of bosons and fermions are fixed) which may 
be identified with the Bosc and Fermi chemical potentials 
and nl and are the values of Ub and n/ at the min- 
imum. Setting the derivatives of £ with respect to 
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and TT-b equal to zero, one finds: 



dn 



f 



\2n 



gbfUf , 
(15) 



(16) 

The solution of Eqs. (fT6|) and p5|) . which are exactly 
Eqs. PT|) and if one drops out the space derivatives, 
gives a minimum if the corresponding Hessian of £ is 
positive, i.e. if: 



dnl duj: \dnbdnf J 



>0, 



(17) 



The solution of this inequality gives the region in the 
parameters' space where the homogeneous mixed phase 
is energetically stable. 

In Fig. [l]we show the region of stability of the homo- 
geneous mixture by the shaded (gray) area in the plane 
{ni,,nf) for different values of ghf. Note that the sign of 
gi,f is not important because in Eq. (jl7p appears g^j:. 
The figure shows that, by increasing g^f, the uniform 
mixture is stable only at large fermionic densities Uf. 
This result, that is in agreement with previous ID pre- 
dictions [3J| on a ID mixture of Bose-condensed atoms 
and normal Fermi atoms, is exactly the opposite of what 
happens in a 3D mixture of bosons and fermions. In fact, 
a uniform 3D mixture is stable only for small values of 
the fermionic density [H, HO, [Ml, [sl, IH, [H]. Figure 1 
also shows that when bosons enter in the TG regime (i.e. 
gb/nb 3> 1) the uniform mixture is much less stable: see 
the behavior of the stability line at very small rib, with 
9bf^0. 

The panels of Fig. 1 suggest that, with a finite gbf, at 
small fermionic densities the uniform mixture is unsta- 
ble: the ground-state of the system displays demixing if 
gbf > and becomes a localized Bose- Fermi bright soli- 
ton if (?b/ < 0. These effects are clearly shown in Fig. [2l 
where we plot the density profiles nj{z), with j = b,f, 
for different values of gbf calculated by directly solving 
Eqs. (jlip and ^12\i numerically with appropriate bound- 
ary conditions. 

The profiles of probability densities plotted in Fig. [2] 
demonstrate the uniform-to-localizcd and the uniform- 
to-dcmixcd transition. In fact, a large and negative 
gbf induces a strong localization while a large and posi- 
tive gbf produces demixing. These density profiles have 
been obtained by solving Eqs. ^ and (fTU|) with a 
finite-difference Crank-Nicholson algorithm and imagi- 
nary time [s^l . 

In the regime gb/rib ^ 1 (BEC limit) and \gf\/nf <C 1 
(BCS limit), the above analysis yields simple analytic 
results. In this case £ = gbnl/2 + 7r^nj/12 -I- gbfnbUf, 
consequently, Eq. ([T7]) leads to the condition 
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FIG. 2: Probability density profiles ni,{z) of bosons (dashed 
line) and n / (z) of fermions (solid line) of the Bose- Fermi mix- 
ture with Nt — 100 bosons and Nf — 10 fermions and peri- 
odic boundary conditions (axial length L = 2ttR = 10). We 
choose = 0.01 and (?/ = —0.025 and calculate the profiles 
for different values of the Bose- Fermi interaction strength gt / . 



for the stability of the uniform mixture. The condi- 
tion Ti^gbnf/2 > g^j is qualitatively consistent with 
Fig. [T] (a). For example, for gb = 0.1 and 5/ —0.1 
the condition of stability of uniform mixture becomes 
Uf > 2glf and for \gbf \ — 0,0.2,0.4, and 0.6 leads to 
iif > 0, 0.08, 0.32, and 0.72 respectively. 

In the regime gb/rib -^i 1 (BEC limit) and \gf/nf\ ^ 
1 (unitarity or tightly bound molecule formation limit) 
£ — gbnl/2 + n^n^/48 + gbfUbUf, consequently, Eq. (fTT)) 
leads to the condition 



-TT gbUf > g^j. 



(19) 



for the stability of the uniform mixture. For exam- 
ple, for gb = 0.1 and gf = —30 the condition of sta- 
bility of uniform mixture becomes Uf > Sg^j: and for 
1.96/1 = 0,0.1,0.2, and 0.3 leads to Uf > 0,0.08,0.32, and 
0.72, respectively, in qualitative agreement with Fig. [T] 
(b). The time-independent conditions ^5]) and are 
also derived in following subsection using rigorous time- 
dependent modulational instability analysis of the uni- 
form mixture (constant-amplitude solutions) in the weak 
and strong Fermi- Fermi coupling limits. 

Finally, in the regime gb/nb 3> 1 (TG limit) and 
\gf/nf\ < 1 (BCS limit) £ = 7r^nl/3 + TT^n)/l2 + 
gbfUbUf, consequently, Eq. pT]) leads to the condition 



TT^nbnf > glj. 



(20) 



for the stability of the uniform mixture. The functional 
dependence of condition ([20]) on nb,nf,gb and g/ is qual- 
itatively different from conditions (US]) and pT ^ 
The inequality (|17[) can be written as 



CbCf > 



IbUf 



(21) 
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where c/ = y'2n/(9/x//9n/) is the sound velocity of the 

superfluid Fermi component and Cb = 2rn,(d / dru,) is 
the sound velocity of the Bose gas. The sound velocity 
Cbf of the ID Fermi-Bosc mixture can be easily obtained 
following a procedure similar to the one suggested by 
Alexandrov and Kabanov [s^ for a two-component BEC. 
One finds: 



Cbf 



1 

71 



c)Y + lQgl^nbnf . (22) 



Thus the sound velocity has two branches and the ho- 
mogeneous mixture becomes dynamically unstable when 
the lower branch becomes imaginary. 



B. Modulational Instability 

We show in the following that for attractive Bose- Fermi 
interaction, the transition from uniform mixture to local- 
ized soliton-like structures considered in Sec. IIIA is due 
to modulational instability To study analytically the 
modulational instability [53, [gO] of Eqs. (O and pU]) we 
consider the special case of weak Bose-Bose (small posi- 
tive gb) and both weak (small <SC 1 corresponding 
to BCS hmit) and strong (large ^ 1 correspond- 
ing to molecular unitarity limit) Fermi- Fermi interactions 
while these equations reduce to 



(23) 



idti^f = [-dl + - gbim""] ^} , (24) 

where k = 1/4 for \gf/nf\ <C 1 and 1/12 for \gf/nf\ ^ 
1. Here we have taken the interspecies interaction to be 
attractive by inserting an explicit negative sign in gbf. 

We analyze the modulational instability of a constant- 
amplitude solution corresponding to a uniform mixture 
in coupled equations and ([M)) by considering the 
solutions 



^60 = AboeM^Sb) = Aboe^'^s^i'^^f^-s^'^io)^ 



Vfo = Afo exp(i(5/) 



Afoe'*'^s-f^'^o- 



(25) 
(26) 



of Eqs. ([22]) and respectively, where Ajq is the 

amplitude and Sj a phase for component j = b, f. The 
constant-amplitude solution develops an amplitude de- 
pendent phase on time evolution. We consider a small 
perturbation Aj cxp{iSj) to these solutions via 



ifij = {Ajo+Aj)c^p{iSj), 



(27) 



where Aj = Aj{z,t). Substituting these perturbed solu- 
tions in Eqs. and and for small perturbations 
retaining only the linear terms in Aj we get 

i dtAb + dlAb - gbAlo {Ab + Al ) 

+ gbfAboAfoiAf +A}) ^0, (28) 

i dtAf + dlAf - 2KT^^A)^{Af + A)) 

+ gbfAbQAfo{Af+A})^Q. (29) 



We consider the complex plane-wave perturbation 

Aj {z,t) ^ Aji cos{Kt - rtz) + iAj2 sm{Kt - nz) (30) 

with j = 6, /, where Aji and Aj2 are the amplitudes for 
the real and imaginary parts, respectively, and K and f2 
are frequency and wave numbers. 

Substituting Eq. ((SQl) in Eqs. and and sepa- 
rating the real and imaginary parts we get 

-AbiK = Ab2n^, (31) 
-Ab2K = Abi^^ - 2gbfAboAfoAfi + 2gbAloAbi, 

(32) 

for j = h, and 

-AfiK = ^/2^^^ (33) 
-Af2K = AfiVl?- -2gbfAboAfoAbi+^KTi'^A'^fQAfi, 

(34) 

for j = f. Eliminating Ab2 between Eqs. (fSTI) and ([32|l 

we get 

Abi[K^ - n^n^ + 2gbAlo)] = -2AfigbfAboAfon^,i35) 

and eliminating Af2 between (|33p and ((34)) we have 

Afi[K^ - + iKTT^A%)] = -2AbigbfAboAfon%36) 

Finally, eliminating Abi and Afi from (j35p and ()36p . we 
obtain the following dispersion relation 



K 



(O^ + gbAl„ + 2ti7r^A%) 



± 



{{9b A 



2KTT^Aj„) 



^9bf^bO^}o 



1/2 



1/2 



(37) 



For stability of the plane-wave perturbation, K has to 
be real. For any Q this happens for 

v2 



igbAlo + 2,i7r^Ajof > (gbA 



\2 

Ho 



2k'k''^ A^-Q j 



^9bfAboA^fQ, 



(38) 



or for 



2gbKT:^Ajf) > glf. 



(39) 

However, for 2gbK'!T^ < gf^^, K can become imaginary 
and the plane- wave perturbation can grow exponentially 
with time. This is the domain of modulational insta- 
bility of a constant-amplitude solution (uniform mixture 
of Sec. HIA) signaUing the possibility of coupled Bose- 
Fermi bright soliton to appear. Noting that = nj 
for the uniform mixture, this result is consistent with 
the stability analysis of Sec. HIA. In the weak-coupling 
BCS hmit k = 1/4 and Eq. ^ reduces to Eq. ([THD 
obtained from energetic considerations. In the strong 
coupling unitarity limit k = 1/12 and Eq. (j39p reduces 
to Eq. (fro)) for the stability of the uniform mixture. Fi- 
nally, we comment that Eq. (j^O)) can also be derived in a 
straightforward fashion from the modulational instabil- 
ity analysis of an appropriate set of dynamical equations, 
with (7b I -06 P replaced by 7r^|-0f,|'^ in Eq- (ESj) and k = 1/4 
in Eq. 



7 



IV. MIXTURE WITH OPEN BOUNDARY 
CONDITIONS 

In this section we consider the ID mixture with open 
boundary conditions. These boundary conditions can be 
used to model a quasi-lD cyHnder. In practice, if the ax- 
ial length L of the atomic wave-guide is much larger than 
the transverse width, i.e. L 3> a± j, j ~ b, /, then one has 
a quasi-lD cylinder. The quasi-lD cylinder becomes in- 
finite for L — > oo. It this case we can use the Lagrangian 
density ([T]) with z £ {—oo, -|-oo). With a finite number of 
particles (bosons and fermions) the uniform mixture in a 
infinite cylinder cannot exist, but localized solutions are 
instead possible with an attractive Bose-Fermi strength 
i9bf <0)M- 



the term involving the LL G{x) function can be evalu- 
ated analytically in the GP {x 0) and the TG (a; ^ 1) 
limits. The analytic term involving the G{x) function 
in Eq. (j43p is exact in both the limits and provides a 
good description of the integral at other values of argu- 
ment provided we take for the constant c = ■\/37r/2. The 
fermionic integral over the GY function F{x) is evalu- 
ated similarly. There is no obvious reason for choosing 
the constant d in this integral. A fter a little experimen- 
tation it is taken to be d = c = •y/37r/2 which provides a 
faithful analytical representation of the numerical result. 

The first variational equations emerging from Eq. ((43|) 
dL/dfit = dL/dfif = yield At = Aj ^ 1. Therefore 
the conditions A^ ^ Af ^ 1 will be substituted in the 
subsequent variational equations. The variational equa- 
tions dL/dwj = lead to 



Variational Results 



Here we develop a variational localized solution to Eqs. 
(fTTj) and (fT2|) noting that these equations can be derived 
from the Lagrangian [6l| 



•/ —OO L 

- <P]F{gf/2<p)) - g,f^l4>) 



dz- ^lbNb- fXfNf (40) 



V 2A^b / 2V2^ V 2Nb J 



dbjNfwl 



^2)3/2 



= 0, 



I ^^/p ( cgfWf \ _ NfgfWf , / cgfWf \ 



gbfNbwj 



0. 



(44) 



(45) 



by demanding 6L/S(j)b = 5L/54>f = 5L/S^,b = SL/Sfif = 
0. 

To develop the variational approximation we use the 
following Gaussian ansatz [6^ 



- ^-1/4, 



-1/4, 



iNbAb 

Wb 

Wf 



■ exp 



exp 



2^L 



2Z.2 



(41) 
(42) 



where the variational parameters are Aj, the solitons' 
norm, and Wj widths, in addition to fij. Note that this 
Gaussian ansatz can be extended to include time de- 
pendence [11], as done recently to investigate the collec- 
tive oscillation of a quasi-lD mixture made of condensed 
bosons and normal fermions (5^ . 

The substitution of this variational ansatz in La- 
grangian (1^0)) yields 



L = fibNbAb + HfNfAf 



NbAb NfAf 



2wl 



2^) 



•wl V 2NbAb J a/Sttw^. 



gbfNbNfAbAf 



^ /^/ p ( dgfWf \ 
- 2 \2NfAfJ 



(43) 



with c = ^37r/2. The integrals over the G{x) and F{x) 
functions cannot be evaluated exactly for all x. However, 



The remaining variational equations are dL/dAj 
which yield as a function of Wj's, and g's: 



= 0, 



Mb 



1 



Nbgb 
2wb^ 



, ( cgbWb \ 
27r I 2Nb J 



V3Nl f cgbWb \ 



gbfNf 



(46) 



2wj 



TTW. 



/ V 2iv, ; 

Nf9f p, f cgfWf \ 
2wfV2^ V 2A^/ / 



9bfNb 



■ (47) 



Equations p4)) — p7)) are the variational results which 
we shall use in our study of bright Bose-Fermi solitons. 



B. Numerical Results 

For stationary solutions we solve time-independent 
Eqs. (jlip and (|12[) by using an imaginary time prop- 
agation method based on the finite-difference Crank- 
Nicholson discretization scheme of time-dependent Eqs. 
([5]) and (jlOp . The non-equilibrium dynamics from an 
initial stationary state is studied by solving the time- 
dependent Eqs. (O and (dU]) with real time propagation 
by using as initial input the solution obtained by the 
imaginary time propagation method. The reason for this 
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FIG. 3; (Color online) Probability densities from numeri- 
cal solution of Eqs. (|9]) and (|10p (here normalized to unity: 
J_^ni{z)dz — 1) compared with variational results given 
by Eqs. dH]) and ^ for iV^ = 300, .g;, = 0.01,5b/ = -0.4 
and (a) Nf = 20,5/ = -0-1, (b) Nf = 15,5/ = -1, (c) 
Nf = 20,5/ = -10, and (d) Nf = 100,5/ = -10- Of these, 
(a) represents fermions in BCS regime, (b) represents fermions 
in the BCS-to-unitarity crossover, (c) represents fermions in 
unitarity regime, and (d) represents bosons and fermions in a 
partially demixed configuration. 



mixed treatment is that the imaginary time propagation 
method deals with real variables only and provides very 
accurate solution of the stationary problem at low com- 
putational cost [131 ■ In the finite-difference discretization 
we use space step of 0.025 and time step of 0.0005. 
First we report results for stationary profiles of the lo- 



0.28 
0.24 

0.2 
0.16 
0.12 
0.08 
0.04 
0. 

(a) 





FIG. 4: (Color online) Dynamics of the probability density 
profiles of (a) Bose and (b) Fermi solitons of Fig |3] (a) when 
at t = 20 they are subject to a quite strong perturbation by 
setting (f)j{z,t) = 1.1 x (f)j{z,t). The solitons undergo stable 
propagation as long as we could continue numerical simula- 
tion. The initial soliton profile is calculated with imaginary 
time propagation algorithm and the dynamics studied with 
real time propagation algorithm. The soliton profiles are nor- 
malized to unity: nj{z,t)dz = 1. 



calized Bose and Fermi solitons formed in the presence of 
attractive Bose-Fermi and Fermi-Fermi interactions and 
repulsive Bose-Bose interactions. In the presence of weak 
attractive Fermi-Fermi interactions, the fermions form a 
BCS state which satisfy a nonlinear Schrodinger equa- 
tion with repulsive (self defocussing) quintic nonlinearity. 
As the strength of the attractive Fermi-Fcrmi interaction 
increases the fermions pass from the BCS regime to the 
unitarity regime which is described by another nonlin- 
ear Schrodinger equation with repulsive (self defocussing) 
nonlinearity. Hence fermions cannot form a bright soli- 
ton by itself. However, they can form a bright soliton in 
the presence of an attractive Bose-Fermi interaction. 

If the fermionic repulsive nonlinearity is not very large, 
bosons and fermions form mostly overlapping (mixed) 
solitons both in the BCS and unitarity regime. Note 
that in the BCS and unitarity regime the fermionic sys- 
tem becomes repulsive in the presence of Fermi-Fermi 
attraction. However, as the fermionic repulsive nonlin- 
earity turns large the fermionic profile comes out of the 
bosonic profile and partially demixed solitons are created. 
We studied the numerically calculated soliton profiles for 
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FIG. 5: (Color online) Chemical potentials for (a) bosons 
and (b) fermions vs. Fermi- Fermi interaction strength gf for 
Nf = 20, Nb = 300, Qb = 0.01 for g^f = -0.4 and -0.8 
calculated numerically (labelled "num" shown by symbols) 
compared with variational results (labelled "var" shown by 
continuous lines). 



wide range of Bose-Bose, Bose-Fermi, and Fermi-Fermi 
interactions and boson and fermion numbers solving Eqs. 

and (|12p by the technique of imaginary time propa- 
gation and compared them with the Gaussian variational 
results obtained from Eqs. and Except in 

the case of strong demixing, when the fermion profile 
strongly deviates from the Gaussian shape, the agree- 
ment between variational and numerical profiles is quite 
good. 

In Fig. [3] we present typical soliton profiles illustrating 
the change in the results during BCS-unitarity crossover 
as well as the demixed profiles. In Fig. [3] (a), (b), and 
(c) we show the soliton profiles for weak (BCS phase), 
moderate, and strong (unitarity regime) Fermi-Fermi at- 
traction corresponding to strengths gf = -0.1,-1 and 
-10, respectively, for gi, = 0.01, gbf = -0.4, iV/ = 20, 
and Ni, = 300 and compare these with the variational 
results. Figure [3] (d) illustrates a demixed state obtained 
by increasing the fermion number from the configuration 
of Fig. [3] (c) from Nf = 20 to 100. In this case the 
(numerically obtained) fermion profile stretches far be- 
yond the bosonic profile and is poorly represented by a 
Gaussian shape, which is the cause of deviation of the 
variational result from the numerical result. 

After illustrating the soliton profiles in different states 
it is now pertinent to verify if these solitons are dynami- 
cally stable under perturbation. To this end we consider 
the typical stationary soliton of Fig. [3] (a) (obtained by 
the imaginary time propagation method) and inflict on it 



the perturbation 4>j{z, t) = 1.1 x 4'j{z, t) and observe the 
resultant dynamics (obtained by the real time propaga- 
tion method), which is illustrated in Fig. [H The solitons 
under this perturbation execute breathing oscillation and 
propagate for as long as the numerical simulation was 
continued without being destroyed. This demonstrates 
the stability of the solitons under perturbation. 

Finally, in Figs. O we show the chemical potential for 
/ib bosons and / for fermions as function of Fermi- Fermi 
interaction strength gf for Nf = 20, N^ = 300, gb = 0.01, 
and gbf = —.4 and —0.8 obtained from numerical solu- 
tion and Gaussian variational analysis. The agreement 
between the two results for ^b is good whereas for jjif 
is only fair. The increase of the Fermi- Fermi attraction 
strength \gf\ for a fixed gbf corresponds to a reduction in 
both chemical potentials signaling strongly bound soli- 
tons. The same happens for the increase of Bose-Fermi 
attraction strength \gbf \ from 0.4 to 0.8. In Figs. [5] the 
small \gf \ limit corresponds to the BCS phase of bosons 
whereas the large \gf\ limit corresponds to the unitar- 
ity regime of fermions. The intermediate values of \gf\ 
denote the crossover from BCS to the unitarity regime. 



V. SUMMARY AND CONCLUSION 

In this paper we have studied a one-dimensional super- 
fiuid Bose-Fermi mixture using a set of coupled mean- 
field equations derived as the Euler-Lagrange equation 
employing the Lieb-Liniger energy density of bosons and 
the Gaudin-Yang energy density of fermions and an in- 
teraction term proportional to the product of probability 
density of bosons and fermions. This set of coupled equa- 
tions has a complex nonlinearity structure for fermions 
and bosons and shows the proper transition from a cubic 
Bose nonlinearity in the weak-coupling GF BCS limit to 
a quintic nonlinearity in the strong-coupling TG (Tonks- 
Girardeau) limit. In addition, for fermions with attrac- 
tive interactions considered in this paper, it shows the 
proper transition from the weak-coupling BCS regime to 
unitarity limit: both limits are described by a quintic 
nonlinearity with different coefficients. In this model, 
in the extreme weak-coupling BCS limit the superfluid 
fermionic energy density is identical to that of a non in- 
teracting degenerate Fermi gas in the normal state [s^ . 

We consider two distinct situations for the superfluid 
Bose-Fermi mixture: (i) in a ring with periodic bound- 
ary condition realizable from a toroidal trap in the limit 
of strong transverse conflnement and (ii) in an inflnite 
cylinder with open boundary condition realizable from 
an axially-symmetric trap in the limit of strong trans- 
verse and zero axial conflnement. 

For the mixture in a ring, from energetic considera- 
tions, we obtain the condition of stability of a uniform 
mixture with a constant probability density. The uniform 
mixture is energetically stable for interspecies attraction 
strength \gbf\ below a critical value, above which stable 
lowest-energy states are bright Bose-Fermi solitons. For 
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repulsive interspecies interaction the stable lowest-energy 
states are demixed states of Bosc-Fcrmi mixture, where 
the region of a maximum of boson probability density 
corresponds to a minimum of fermion probability den- 
sity. It is also demonstrated algebraically that for at- 
tractive Bose-Fermi interaction the bright solitons can 
be created from the uniform mixture above the critical 
Bose-Fermi attraction by modulational instability of the 
uniform mixture under a weak perturbation. 

In the one-dimensional infinite cylinder we solved the 
coupled set of equations for the superfluid Bose-Fermi 
mixture numerically and using a Gaussian variational 
approximation. We calculated numerically the probabil- 
ity density profiles of the bright solitons as well as their 
chemical potentials and compared them with the respec- 
tive Gaussian variational approximations. The agree- 
ment between the two is good to fair. In this case a 
partial demixing of the Bose-Fermi solitons is possible, 
when the Fermi soliton extends over a very large region 
in space while the Bose soliton remains fairly localized. 



We also established numerically the dynamical stability 
of the Bose-Fermi solitons by infiicting a perturbation 
on the solitons by multiplying the wave-function profiles 
by 1.1. The system is then found to propagate over a 
very long period of time executing breathing oscillation 
without being destroyed, which demonstrated the stabil- 
ity of the solitons. Finally, we comment that in view 
of the experimental realization of the superfluid Bose- 
Fermi mixture psj and observation of solitons in a pure 
BEG [36l . [13, the achievement of a bright Bose-Fermi 
soliton seems possible through a controlled manipulation 
of strengths of atomic interactions by varying an exter- 
nal background magnetic field near a Feshbach resonance 
and by adopting the set of parameters we have used 
in this paper. 
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the Italian GNFM-INdAM through the project "Giovani 
Ricercatori" and by Fondazionc GARIPARO. 
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